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Abstract. The transition from acoustic noise in the radiation-dominated universe to the density structures in the 
matter dominated epoch is considered. The initial state is a stochastic field of sound waves moving in different 
directions. The construction of the initial state is compatible with the hyperbolic type of propagation equation for 
density perturbations, and parallel to the theory of stochastic background of gravitational waves. Instantaneous 
transition between the cosmological epochs is assumed, and Darmois-Israel joining conditions are applied to match 
solutions for sound waves with growing or decaying modes at the decoupling. As a result a substantial amplification 
of the low scale structures is obtained. 
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1. Introduction 

Simple models of transitions between different cos- 
mological epochs help us to understand amplifica- 
tion of scalar fields, electromagnetic and gravitational 
waves in the expanding universe (Frieman & Turner 
1984, lAbbott fc Harari 19861 IHu 19981 ICrishchuk 19961 
I Allen et al. 20 00 ) . Typical cosmology of that class consist 
of three phases: 1) a semiclassical phase, which is com- 
monly identified with de Sitter stage, 2) the radiation- 
dominated epoch described by the equation of state p = 
e/3, and 3) the matter dominated era when pressure is 
negligible p — 0. Transitions between them are assumed 
to be instantaneous. The physical meaning of these models 
is close to that of simple models of the particle scattering 
on rectangular potential barriers in quantum mechanics. 
Asymptotic results weakly depend on the barrier profile, 
so we hope, that in more realistic cosmology, the pertur- 
bation amplitude in the remote past and the far future 
should only marginally depend on the transition details. 

For scalar or electrodynamic fields, as well as for 
gravitational waves the amplification can be measured 
by Bogolubov coefficients (IBirrell fc Davies 1982J) . These 
fields arc governed by hyperbolic partial differential equa- 
tions in each of discussed epochs and the changes in the 
equation of state - the change in the background dynam- 



ics - result in a specific relations between Fourier modes 
in the in and out state, respectively. For some frequen- 
cies, counter-propagating waves are substantially ampli- 
fied (creation of pairs of particles with opposite momen- 
tum) - the so called back scattering effect (Parker 1972) 
comes into play 1 . Similar effects may occur in acoustic 
field IjLukash 1999|l . Fields enhancement on non static 
background is generally referred to parametric amplifica- 
tion HGrishchuk 1995jl . 

Neither parametric amplification nor particle cre- 
ation theories provide appropriate language to describe 
the growth of density perturbations in the transi- 
tion to matter-dominated epoch. The scalar pertur- 
bations form waves llSachs fc Wolfe 19671 IWhite 19731 
Field & Shepley 1968) (in quantum phonon-approach 
iLukash 19801 IChibisov fc Mukhanov 1982|> only in the 
epochs of non vanishing pressure, but they transform 
into non-travelling ingomogencities (growing and decay- 
ing modes) in the matter dominated epoch (p = 0). 
The change then is two-fold: 1) the transition modifies 
the background dynamics, 2) the propagation equation 
change its differential type. Bogolubov coefficients loose 
their physical meaning. Yet, the general scheme of the 
field propagation throughout the transition epoch is the 
same: for the scalar or electrodynamic fields the continuity 
of each field and its time derivative must be satisfied at 
the transition, for metric perturbations (both scalar and 



1 For low frequency gravitational waves standing waves of 
substantial amplitude may appear (squeezed state of gravita- 



tional field (Gr ishcliuk 
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tensor) the first and the second fundamental forms must 
be continuous (Darmois-Israel conditions). 

Classical perturbations in the cosmological mod- 
els with sharp transitions have been investigated by 
Kodama and Sasaki (Kodama & Sasak i 1984|) or Hwang 
and Vishniac | |Hwang fe Vishniac 199l| . These authors 
limit themselves to the standing wave solutions and to 
regime of low frequencies. A similar task for a uni- 
verse with the radiation and dust mixture has been 
undertaken by Mukhanov, Feldman and Brandcnberger 
(Mukhano v et al. 1992f) also in the long wave limit. In 
quantum theories the radiation to matter transition has 
been extensively discussed by Grishchuk ( Grishchuk 1994) 
as a part of a more complex cosmological model including 
an inflationary stage. Large scale perturbations where in- 
vestigated there in the context CMBR temperature fluc- 
tuations. Relatively little is known about the low scale 
inhomogeneities . 

In this paper we present exact formulae for density 
perturbations in the universe with sharp transition in a 
full range of frequencies. We start with acoustic noise in 
the radiation era and investigate its transition to density 
structures in the matter dominated epoch. We employ 
the autocorrelation function as a measure of structure, or 
equivalently its Fourier transform-the spatial spectrum 
of inhomogeneities. This measure agrees well with what 
other authors propose in cosmology. However, despite the 
cosmological practice we do not limit the basis of elemen- 
tary solutions to " growing modes" . Instead, we take into 
account the complete basis of Fourier modes, adequate 
to the hyperbolic character of the propagation equation 
and in full analogy to a stochastic description of gravita- 
tional waves IjAllen et al. 2000)) and quantum (phonon) 
theories l|Lukash 19801 IChibisov fc Mukhanov 19821 
IGrishchuk 1991 . Since we disregard inflation, there 
is no squeezing mechanism (Grishchuk 1995) for the 
acoustic field. We examine the spatial power spectrum 
without being limited to standing waves, and show that 
substantial amplification occurs in the high frequency 
regime, which has not been investigated till now. 

The structure of the paper is the following: in Section|21 
we discuss exact solutions to the perturbation equations 
in both radiation and matter eras employing orthogonal 
gauge. In the Section we construct the 2-epoch model 
with the instantaneous transition in the equation of state, 
by use of the Darmois-Israel conditions. The Section 0] is 
devoted to the stochastic properties of the acoustic field. 
Finally, we investigate the time evolution of the spatial 
spectrum of cosmological inhomogeneities. 



2. Perturbations in the Robertson-Walker universe 

There are several independent methods to de- 
scribe scalar perturbations in a gauge-invariant way 
l|01son 19761 IBardeenT980l [ Brandenberger at al. 1983| 
|Lyth fe Stewart 19901 lEllis fc Bruni 19890 . For the uni- 



verse filled with an ideal fluid 2 with arbitrary non- 
vanishing pressure p all of them lead to the same 
result: under appropriate choice of the perturbation 
variables, the propagation equations converge to the wave 
equation (|Sa,chs fc Wolfe 19671 |Field fc Shepley 19681 
IChibisov fc Mukhanov 19821 |Golda fc Woszczyna 200iy 
The density perturbations form travelling waves. Below 
we consider rotationless fluid, hence the hypersurfaces 
orthogonal to the fluid flow can be globally defined, and 
therefore, the orthogo nal gauge (|Lyth fc Mukherjee 1 988 
|Lyth fc Stewart 1990| IPad manabhan 1993) is naturally 
applied. 

The linear corrections to the energy density and the 
expansion rate evolve according to 

d t Se(t, x) = -p 5&{t, x) - e 5ti{t, x) - # 5e(t, x) (1) 

dt60(t, x) = -4GirSe{t, x) - V ^' x) (2) 

Po + £o 

- Wx)- (24G7r r w,x) 

3 3(p + eo) 

where the subscript "0" refers to the unperturbed, back- 
ground FRW universe (eo = eo(i), and po = po{t) are 
solely functions of time). The equations £Q), (0) are de- 
rived from the Raychaudhuri and the continuity equations 
by use of the standard linearization procedure and by re- 
placing the proper time by the orthogonal time. Below we 
limit ourselves to spatially flat universe K — 0. In this case 
the background energy density and the expansion rate are 
related to each other by 87rGeo = $ 2 /3. 

Let us assume now that the universe is filled by the sin- 
gle fluid with the equation of state P/e = w = constant, 
where the value of w determines both the evolution of the 
background metric (unperturbed universe) and the sound 
velocity for small perturbations. We express the energy 
density and the expansion rate as the composition of the 
background energy values eo, $o and the small inhomoge- 
neous correction Se — eo^(i,x), S$ = $o©(^ x ) 

e(i,x) = e (t)(l + <5(t,x)) (3) 

0(t,x) = o (t)(l + e(t,x)) (4) 

where S(t,x) and 9(i,x) play the role of the density 
and expansion contrasts, respectively. Transforming sys- 
tem (j2J) to a second order propagation equation for 8 
we obtain a partial differential equation of the form 

wL35(t,x) = -?-(w- l)(l + 3wWi,x) 
6 

+ 4^-»)w,x)+3 f 2 J(i,x) (5) 

with the background evolution given by $o (t) — (i+ w )t or 
equivalently in the conformal time r\ 

wL3 ^' x) = (£tt)^ (7? ' x) 

- (^^j-dAv^ + d^x) (6) 

\ 6w + 1 / 7] ' 

2 With the diagonal energy momentum tensor T M „ = (e + 
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2uj \ lUT) 



= i(kx-wij) 



(12) 



The conformal time rj is defined here as the integral The modes u k ( e )(7j,x) are simply -p==e 4 ( kx Ljr '\ while 
/ W) dt of tne scale factor reciprocal over the orthogo- Uk(e) (jj !X ) can be found as 
nal time t and L3 stands for the Laplace operator in 3- 
dimensional Euclidean space. Equation JHJ can be solved 
analytically by use of Fourier transform. Actually, we are 
interested in two special cases w — 1/3 (radiation- filled 
universe) and w = (matter-domination epoch). 

The propagation equations for small perturbations , 
(0 (and consequently Q below) do not contain the gravi- 
tational constant G, which means that inhomogeneities do 
not self-gravitate unless the linear regime breaks down. 
They evolve as acoustic waves in the expanding gas 
medium (compare ISachs fc Wolfe 19671 IStoile 2000j> V All 
the perturbation equations obtained in different gauge in- 
variant formalisms can be reduced to an equation of the 
form JSJ by suitable changes of variables (and with dif- 
ferent meaning of the variable S and the time parame- 
ter) ( |Golda fc Woszczyna 200T] ). The necessary changes 
reads (in the notation: reference: original notation — > 8) 
(ISakai 196911: K -> 6; 
HBardeen 1980J) 



5; 



Pm — > S] 

IjKodama fe Sasaki 1981 c hap. IV): A 
I Lyth fe Mukherjee 1988| ): S -> 6; 
l|Padmanabhan 1993|l : S — > 5; 
flBrandenberger at al. 1983| |: <!>hI pa 2 — > 5; 
|Ellis et al. 1990jl : V ->lL 

Transformations of these equations to conformal time (if 
parameterized differently) are necessary. 

In the universe filled with highly relativistic matter, 
the scale factor a(rf) evolves as a linear function of the 
conformal time: a(rj) = y^M/3r], and preserve M. = e^a 
as the constant of motion. Equation JfJJl expressed in con- 
formal time takes the canonical form (independent of first 
derivatives) 



The generic perturbation <5(?7,x) is composed of 
travelling plane waves u k ( e j with decreasing ampli- 
tude. (Similar solutions are known in the theory 
of gravitational waves IjWhite 19 92). and scalar field 
|Stebbins fe Veerarghavan 1993| ). The Fourier coefficient 
A k = — iujAk is an arbitrary complex function of the wave 
number k, while 5(ry,x) keeps real values. The frequency 
uj obeys the dispersion relation uj 2 = k 2 /3, hence waves 
of all length-scales propagate with the same phase and 
group velocity 3 (compare IChibisov fc Mukhanov 1 982 
ISachs fc Wolfe 19671 IWhite 1973|> . Modes u k(e) , although 
different from the simple eikonal form, are still orthonor- 
mal in the sense of the Klein-Gordon scalar product. 

In the epoch of matter dominance m = cqo 3 is the 
constant of motion and the scale factor evolves as a(rj) = 

2 

^jo - • The propagation equation expressed in orthogonal 
gauge (in all formalism mentioned above) reads 



G 
rj 



<5(r?,x) + -d v 6(r), x) + d 2 6(r/,x) = 0. 



V 



(13) 



1 r L35( ?7 ,x) = -%^ + ^( ?7 ,x) 



(7) 



3 rr 
Equation J7J is essentially the same as the propagation 
equation for gravitational waves in the dust-filled uni- 
verse (|Grishchuk 19741 IWhite 1992|) . It reduces to the 
wave equation in its normal form 

\L35(r,,x) = dfi( V ,x.) (8) 



for S(n,x) defined as 
%, x ) = ~d n (r)5(r),x)), 

7] 



(9) 



The variable <5(?7,x) is the orthogonal-gauge analogue to 
the Field- Shepley variable H ( |Field fc Shepley 1968|), 
or the Sachs-Wolfe variable E ISachs fc Wolfe 19671 
Solutions S(i],x) and S(n,x) expand into Fourier se- 
ries in the way appropriate for massless scalar fields 
IRirrell fc Oavies 1982|) . 



(10) 



<ty7,x) = / (AkUk( £ )(»?,x) + A£u£ (e) (j7,x))<ik 



<5(?7,x) 



J (4«k( f )(l,x) + .4 k u k(e) (?7,x))dk 



(11) 



Vanishing pressure implies the absence of the Laplace op- 
erator, consequently, the general solution consists of grow- 
ing and decaying solutions involving two arbitrary func- 
tions of the space coordinates fx (x) and f% (x) : 

<%,x) = / 1 (x)7 ? 2 + / 2 (x)rr 3 (14) 
This solution expands into Fourier series 

5(7], x) = J (a k e* kx + a k e~ 4kx )r/ 2 + (fo k e* kx + & k e- ?;kx )?r 3 dk 

(15) 

where the coefficients a k and & k are arbitrary complex 
functions of k. 



3. Matching conditions in the transition epoch 

Consider now the two-epoch cosmological model com- 
posed of both, the radiation epoch (governed by the equa- 
tion of state p = e/3), and the succeeding epoch of matter 
domination (with p — 0). Below, the quantities related to 
these two epochs will appear with the indices (1) and (2), 
respectively. We assume that the transition between the 
epochs is instantaneous and occurs on the hypersurface E 
orthogonal to the four velocity u M of the matter content. 

The initial Cauchy conditions are unique and con- 
sistent on the hypersurface E if the first and the 

3 This is not the case of open or closed universes, 
where sound are dispersed on the space curvature 
HGolda fc Woszczyna 2551) . 
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second fundamental forms are equal IjDarmois 19271 
|Hawking fc Ellis 1973) 



V(i) 



(S) 



V (2) 



(16) 
(17) 



Subscripts (1) and (2) refer to both half-spaces divides by 
the surface £. For unperturbed background these equa- 
tions imply the continuity of the scale factor a(rj) and its 
first derivative 



0(i)(E) = 0( 2 ) (S), 
5^0(i) (S) = 9,,a(2)(S). 



(18) 
(19) 



In our two-epoch model the matching conditions are sat- 
isfied by 



0(1)(7?) = 



4%/3r; s 



(20) 
(21) 



where r/s denotes the time of the transition, and M. is 
the constant of motion M. = £o a (i) of the radiation filled 
universe. 

In the gauge-orthogonal formalism conditions ((161 I17JI 
can be rewritten to directly join the density and expan- 
sion perturbations in both epochs, before and after the 
transition. From 1(161 I17[) one can easily find that the en- 
ergy density e and the expansion rate $ are continuous 
on £. Indeed, the energy density e is related to the in- 
duced curvature R^ and the second fundamental form 
by Pawking fc Ellis 1973) 



2e = R (3) - (x£) 2 - X^t 



The Ricci scalar on £ consists of the metric form 
h^viYi) and its space derivatives, so R^ like /i M „(£), is 
continuous in the transition, Xfiv i s continuous on the 
strength of 117JI . As a consequence, 1)22(1 ensures, that e 
is continuous. On the other hand the continuity of the 
expansion rate $(t) comes directly from the same prop- 
erty of the second fundamental form x M „ . Eventually, af- 
ter eliminating unperturbed values eo and i9 : the density 
perturbations obey 



£tf (1) (£) 
where 
Se(j], x) 
5d{r], x) 



= <fe(2)(£), 
= ^(2)(S) 



Po + eo 



<Je(?7,x) - 



a(P + eo) 



(23) 
(24) 



(25) 

d^efox), (26) 



as derived from and ©. It is important to note that 
the time derivatives of e and $ may jump, therefore neither 
Se'(r]) nor 5'(rj) = (5e(r])/e)' are continuous on £. 

The starting point for the investigations of acoustic 
fields and the structures they produce in the two-epoch 



universe is the analysis of a single monochromatic wave (a 
single Fourier mode) 



«k(e,l) 



1 



1 



ILU7] J 



i(kx — LJT}) 



(27) 



This wave, while falling on E, generates on the "other 
side" of the transition surface a mixture of growing and 
decaying mode 



M k(e,2) — Ok 



,?'kx 



?7s 



(28) 



The two coefficients ak and &k are uniquely determined 
from P3"l 12*3) , Indeed, evaluating J"""3 [2"") for the pertur- 
bation 1(271 I28|l and imposing the matching condition 1(231 
12411 one obtains 



"k = 



40 
1 



12 



(29) 
(30) 



The frequency lo refers to the acoustic wave in the ra- 
diation epoch, and is related to the wave number fc by 
the linear dispersion relation ui — k/y/3. In the epoch of 
matter dominance, the perturbations lose their wave char- 
acter, therefore it may be better to parameterise them by 
the wave number fc, which has a well defined meaning in 
both epochs. Now, the modes Uk(e) with coefficients l|2*9l 
13011 take the form 



Wk(£,l 



(22) «k( e 




(31) 



(32) 



The second independent solution consists of complex con- 
jugates of Wk(e,i) and Wk(e,2)- We will include it later to 
restore the complete Fourier basis. 

4. Random acoustic fields 

The acoustic field in the early universe is shaped by ther- 
modynamic or quantum phenomena acting prior to and 
during the radiation era. Their probabilistic nature leads 
to stochastic description. We restrict ourselves to stochas- 
tic processes homogeneous in the broad sense, called also 
weakly homogeneous processes, which keep their mean 
value and variance (standard deviation) invariant under 
translations. The two-point autocorrelation functions for 
them are functions of the distance between points solely 
- not of these points' positions (Locve 1963). Processes 
homogeneous in the broad sense have their Fourier repre- 
sentations l|Loeve 19631 |Sobczyk 1991[|Yaglom 1961) 



(33) 



8(v,*) = / (AkU k (e)(??,x) + A£u* (e) (?7,x))dk 
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where the integral is understood to be the stochas- 
tic integral, and the Fourier coefficients Ak are random 
variables. Their expectation values fulfil ( |Sobczyk~1 991 
|Yaglom 19611 



E[A k A* k ,} ~ S(k-k') 
E[A k A k ,} = 



(34) 
(35) 



where S denotes Dirac's delta. Conversely, each process 
obeying (|34l I35fl is homogeneous in the broad sense. 

The relations Q34I 13 5|) have a clear physical meaning. 
The first of them expresses the statistical independence of 
waves with different wave- vectors. The second means that 
phases of perturbations at any moment and any place are 
statistically independent. Altogether they assure statisti- 
cal independence of waves moving in different directions. 
This stochastic process can be also expressed in terms of 
the 5(rj, x) variable. With help of the random Fourier co- 
efficients A k satisfying 



E[A k A* k ,} = P k S(k - k') 
E[A k Ai,} = 



we write 



(36) 
(37) 



(38) 



We assume that the power spectrum of the acoustic noise 
Pk depends solely on the magnitude of the wave vector 
and not on its direction. The perturbation is generic - 
no additional constraints, nor any ordering or squeezing 
mechanisms have been introduced. The construction of 
a random acoustic field is identical with the construc- 
tion of random field of gravitational waves in the matter 
dominated universe IjAbbott & Harari 19861 lAllen 1 996 
lAllen fc Romano 19991 |Maggiore 2000| ). It reflects the 
equivalence of the propagation equations for both these 
classes of perturbation. Stochastic acoustic field fulfill- 
ing conditions I|36I37|I can be understood as a classi- 
cal limit of the quantum phonon approach (JLukash 1 980 
isov fc Mukhanov 19821 IGrishchuk 1994|l . 
Modes Uk( e ) while regularly extended across the tran- 
sition hypersurface £ define a stochastic structure (by 
means of integral <|38j) ) in the matter dominated universe. 
The commonly used measure of this structure is the two- 
point autocorrelation function 

R(r), h) = — { E[S(r], xWrj, x + h)Wh-h - 1) dh (39) 
4tt J 

defined on the constant time hypersurfaces. Under condi- 
tions (j34H37(l this function reads 

R(r],h) = -j- I 2u k u k Pfe exp(ik-h)*(h-h - l)dkdh 

47T J 

] - Pk ( V )dk (40) 



, 9 sm(hk) 
47rfc - 



Ilk 



where 



plays the role of the structure spatial spectrum and 
itk( e )(?7,x) are modes defined by l|31l 132(1 extended to both 
epochs before and after the transition. In this way acous- 
tic noise in the radiation era with a given spectrum P k 
determines uniquely (via Darmois-Israel conditions) the 
spatial spectrum of inhomogeneities p k (t)) at any stage of 
the structure formation process. 

It is important to clearly distinguish between p k and 
Pk, and understand their roles in the cosmological con- 
text. As defined above, pk is the Fourier transform of 
the two-point autocorrelation function. In the cosmologi- 
cal literature it is called the power spectrum 4 by analogy 
to similar concept known in the analysis of time series 
l|Anderson 1971|l . but its physical dimension is different. 
On the other hand, P k defined by ffity is the genuine power 
spectrum of the acoustic field (with the same physical 
sense as the Planck power spectrum of the electromagnetic 
radiation), and can be obtained from Hamiltonian descrip- 
tion IjLukash 1980(1 . Although in cosmology one cannot di- 
rectly observe Pk, this quantity defines the physical state 
of acoustic field. This is P k not pk, which should be ei- 
ther guessed, or inferred from fundamental laws of physics 
l|Chibisov fc Mukhanov 1982(1 . The shape of P k is not pre- 
cised in this paper. 



The time factor 



T ( £ ) 



2u 



k( e )W k(e ; 



converts the acous- 



tic spectrum P k into the spatial-spectrum Pk(v)- (It is 
numerically equal to Pkiv) for the white noise acoustic 
field: P k = const.) Factor tia contains the entire time 
dependence of cosmological inhomogeneities. Employing 
extended modes l|3"Tl l3~2*)l . we find 



T (M) 

T (^2) 



t(' + — 

k V k 2 ft 2 
V3/128/ 5 
k I 5 I 2P 



(42) 




3(10- 3A: 2 
25 



-(1 + ^ + 




(43) 



p k {n) = 2M k(e) u£ (e) P fc = T (e) Pk 



(41) 



where fj is the normalized time parameter fj = rj/rfs and 
the k is the modified wave number k = krjs/y/3 which 
measures the number of oscillations within the sound hori- 
zon on S. We hold the factor V3/fc in both Ti e ,i) an d 
r (e,2) to keep modes Uk(e,i) orthonormal in the sense of 
the Klein-Gordon scalar product in the radiation era. It is 
easy to check that the time factor r is continuous on E. 

In the radiation era, perturbations greater than the 
sound horizon decrease, while the low scale ones main- 
tain constant amplitude. This can be seen directly from 
(|42(l . The result qualitatively agrees with quantum the- 
ories fsee IMukhanov et al. 19921 part III, formula 20.6). 
Quantitative difference comes from the different gauge 
choice 5 . It should be emphasized that the spectrum Pk, 

4 In this context we consequently use the name spatial power 
spectrum. 

5 The quadratic behaviour of the decreasing term is char- 
acteristic for the orthogonal gauge. In the synchronous sys- 
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the spatial spectrum pk and the parameter r are invari- 
ant under unitary transformations. Therefore, the result 
is physically well defined - it is not the effect of any par- 
ticular choice of the Fourier basis. 

In the matter era the perturbations (|32|l are a spe- 
cific mixture of growing and decaying solution (compare 
lEllis et al. 19901 |Liang 19771 ). Therefore, despite vanish- 
ing pressure, their evolution depend on their length-scales. 
The dominant growing term (1 + fj) in 1|43(1 is multiplied 
by k 2 , which means strong amplification of short waves. 

A peculiar velocity field can be measured by the de- 
viation 5$ from the homogeneous Hubble flow t9 — 3H. 
Modes of the fluid compression 6 are associated to the den- 
sity modes i[12|). H27(l . On the strength of H26|) we obtain 



50 a 



3^3 



(i) 



2uj 2VMr] 



1 



1 

i-quj 



'_!hL ) e i(kx -.1- , 



(44) 



8{t play a role in the matching conditions, as the Darmois- 
Israel conditions require that both the density and expan- 
sion fields be continuous at the transition. (Their time 
derivatives my undergo discontinuity.) The expansion con- 
trast 9 = S'&f'&o can be expressed as a stochastic integral 

8(77, x) = J (A k u kW (r),x) + Alul m (ri,x))dk (45) 

with the same coefficients Ak as in (|lip. but with a dif- 
ferent modes Uk(#)j given by 

31/4 I 1 , "^3 likr)\ ikx-ifc^= 

u ^ = ^m\ 1+ ^ + -^) e (46) 

3 1/4 lkx -i£2& ( ikriY, A i] \ 2 




(47) 



In close analogy with the density perturbations, one can 
express the peculiar velocities (perturbation in the expan- 
sion rate) in terms of the spatial autocorrelation or spec- 
trum. This time 

R( rh h) = -jjj- J E[G(r], x)0[x + h, r))]5{h-h - 1) dh (48) 

and the time factor r is given by formula (|41|l after re- 
placing u k(e) by u kW 



V3 / fj 2 ~k 2 1 



fj 2 k 2 



(49) 



tem of reference this term decays as - 4 . 4 . It can be eas- 
ily derived in the Field-Shepley formalism, see formula (5.3) 
of ijChibisov fc Mukhanov 19821 after substituting solutions 
of (4.7) and evaluating the integral over 77. Equivalently, 
it can be proved directly in the original Lifshitz formalism 
HGolda fc Woszczyna 2Q0l\ . 

The orthogonal gauge realize the comoving system of ref- 
erence, where the four velocity is globally chosen as u = 
(1,0,0,0), therefore one cannot describe the peculiar velocity 
field directly by Su 




(50) 



Unlike the density spectrum, the spectrum of the expan- 
sion rate monotonically increases in the radiation domi- 
nated universe l|49l see also the behaviour of modes ww^.i) 
— HQ)! . This is so, because the decrease in peculiar veloci- 
ties is less than the decrease in homogeneous Hubble flow. 
Increase of the expansion contrast 8 with time does not 
invalidate the acoustic approximation, as long as the gas 
velocities are small when compare with the sound velocity 
( Whi tham 1974| . In our case the acoustic approximation 
is valid when the contribution of terms quadratic in <5i?(?y) 
to W(i) (0 is negligible. For large values of time this is 
assured by (|44H . Therefore, a perturbation initiated as an 
acoustic perturbation remains acoustic during the entire 
radiation era, even though Stf/fl increases. 

To express the perturbation enhancement relative to 
its initial amplitude we introduce the ratio 



T e (fj,fji,k) 



Pk(Vi) T e (fji,k)' 



(51) 



The so called transfer function 5 /Si is often used in a 
similar context. Fitting to this traditional name we will 
call T(r),k) the spectrum transfer function . Contrary to 
5/5i the spectrum transfer function T is real and invariant 
under unitary transformations of Fourier bases. Plotted 
on a logarithmic scale (Fig. ^ Fig- it shows that the 
perturbations evolution is highly sensitive to k. In both 
Fig. E an d Fig. |21 we start from rji = 0.1 77s (fji = 0.1) to 
evaluate perturbation up to 77 = 3 77s (77 = 3). 

As already mentioned, the large scale inhomogeneities 
(those, which are larger than the local sound horizon k -C 
1) decay in the radiation era (77 < 1). In this case, the term 
■p^- is dominant in the time factor T( e> i) (f), k) and strongly 
decreases with time. Decay of the large scale component 
is a generic feature of acoustic noise on the expanding 
radiation dominated homogeneous background. In other 
words, the homogeneity of the radiation filled universe is 
a stable property at least as long as the generic pertur- 
bations are taken into account. Similar phenomenon of 
large-scale wave extinction is observed in scalar field the- 
ory (Stcbbins & Vccrarghavan 1993). Although the large- 
scale inhomogeneities change substantially during the ra- 
diation era, their response to the change in the equation 
of state is very weak. Their wave character vanishes at 
77s, but the spatial spectrum itself is insensitive to the 
transition (see Fig. [2Jl . 

The low scale perturbations produce a different sce- 
nario. Their amplitudes kept constant during the radi- 
ation dominated-epoch to increase by several orders of 
magnitude at the transition. For inhomogeneities of galac- 
tic scale (M = 10 n M s ) the time factor r is of 10 8 , 
which means a 10 times amplitude enhancement. The 
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log(Te) 




Fig. 1. The spectrum transfer function T e for the density 
perturbations (log scale) as a function of conformal time 
fj. The family of solutions covers the wave numbers range 
log(fc) S [—1,5] i.e. the lowest curve on the diagram refers 
to the perturbation scale ten times larger than the sound 
horizon distance, while the top one to a perturbation scale 
10~ 5 times smaller. Point ij — 1 represents the transition 
from the radiation to the matter dominated universe. 



log(T £ ) 




Fig. 2. The density spectrum transfer function T e as a 
function of both the conformal time i) and the wave num- 
ber k (fj = 1 - the transition from the radiation to the 
matter dominated universe). 

pressure discontinuity excites low-scale perturbations and 
their contribution to the spatial-spectrum becomes domi- 
nant shortly after the transition. 

As opposed to the density perturbations, the veloc- 
ity magnitude does not change significantly at the tran- 
sition. An analogue of the transfer function Tfl(i),k) = 
Tfi{r), fc)/r^(%, k) constructed for the expansion spectrum 
is shown in the Fig. The expansion contrast increases 
systematically during the radiation era and smoothly en- 
ters the matter era, keeping nearly the same growth rate 
before and after the transition. A substantial change in 
the velocity field relates to its phase. Figure Fig. 0] shows 
the difference in arguments of the density and expansion 
modes S(f> ~ arg(itk( e )) — arg(uk(^)). In the radiation era 
large frequency modes, (both Wk( e ) and iik(i9)) are shifted 
to each other by ir/2 (compare also Q46H and 127fl ). This 
effect is the generic perturbation property in the radia- 
tion dominated fluid (Elli s et al. 1990J) . and is an artefact 
of their acoustic character. Similar relation of phases can 
be also obtained on the ground of binding energy analysis 



| |Liang 1977| ). After the transition the arguments of the 
density and velocity modes differ by n so their maxima 
are anti-correlated. Smallest expansion occurs in the re- 
gions of highest density. Decreasing modes have decayed 
and fluid flow becomes potential (jPeebles 1980J) . (Slightly 
more complex is the behaviour of the low frequency modes. 
Solutions like that hardly achieve the phase of potential 
flow (the argument shift is different from 7r).) 



log (To) 




log(k) 



Fig. 3. The expansion spectrum transfer T$ as a function 
of the conformal time 77 and the wave number k. 




log (k) 

Fig. 4. Phases of perturbations. 

^From the hydrodynamic point of view, the transition 
from the radiation to the matter- dominated epoch breaks 
down the acoustic approximation. The sound velocity in- 
stantly falls from v = l/\/3 to zero, and in consequence, 
the fluid velocity formally becomes greater than the sound 
velocity at each point in space. (In this way the struc- 
ture formation has some aspects characteristic for forma- 
tion of acoustic shocks - compare |Montenegro et al. 1999[ 
Kho perskov &: Khrapov 1999| ) The growth of inhomo- 
geneities is based on acoustic instability - the self- 
gravitation processes do not switch on until the linear 
regime fails. In both cases — the density and expan- 
sion fields — the spectrum transfer function monotoni- 
cally increases with the wave number k. In particular 
the phenomenon of acoustic peaks is absent 7 . The ab- 
sence of peaks is characteristic feature of the random 

7 The same effect appears in more sophisticated transition 
models investigated numerically (Press fc Vis hniac 1980), if 
they include the complete basis of solutions. On the other 
hand, as shown in fVoglis 1986), peaks may appear in a 
pure radiation-filled universe model without evoking com- 
plicated recombination processes, if one limits to growing 
modes (standing waves) with specific phase correlation. These 
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field composed of the statistically independent moving 
waves (Grishchuk 1995)) . The spectrum transfer function 
may substantially change in the small fc-regime when the 
universes undergoes more than one phase transition (some 
classical ones, like the change in the sound velocity, or 
semiclassical like the transition from inflationary to ra- 
diation dominated universe). In these cases the squeezed 
states in the acoustic field appear on large scales, which 
are an alternative explanation of the CMBR temperature 
spectrum l)Bose fc Grishchuk 20 02 ). 

5. Growing modes versus complete solutions 

When considering instantaneous transition to the mat- 
ter dominated universe authors limit themselves to the 
large scale perturbation regime. In this regime the grow- 
ing mode amplitude after the phase transition is entirely 
dominated by that before the phase transition and there is 
no chance of generating a growing mode out of a decaying 
mode (|Kodama fc SasakTT 984). Although the short per- 
turbations do not produce the same scenario 8 the decaying 
modes of any length-scale are commonly neglected at the 
beginning of the matter dominated era. 

By abandoning all the decaying terms in the formula 
l|43|l one obtains r e oc (1 + ??) 4 , and consequently, the spa- 
tial spectrum of the density contrast 

Pk(fj) = T e (fj, fji)p k (fji) (52) 

divides into time dependent 

r ^ ) = (rrf) 4 (53) 

and the scale dependent pk{fii) (given at some fji > 1). 
The spectrum transfer function T e drawn in the logarith- 
mic scale is presented in the Fig. GD T e does not depend 
on fc, therefore the perturbation is amplified in the scale- 
independent way. This scale-independence is commonly 
attributed to the dust-filled universe (p=0), but actually 
it requires more restrictive conditions: the perturbation 
must be solely composed of growing modes. 

While the pure growing modes increase ten times in the 
interval ij 6 (1,3), the mixture of growing and decaying 
low scale modes is enhanced by several orders of magni- 
tude (compare Fig. ^ an d Fig. |SJ . The enhancement of 
this mixture depends on the wave number k. Substantial 
magnification of the low scale modes (Fig. ^) means that 
the low-scale structures may enter the nonlinear regime 
first. 

When the decaying modes are taken into account the 
scale- independence breaks down. Decaying modes "re- 
member" perturbations' past. In our case the admixture 
of decaying modes defined by the joining conditions 

phenomena have also been discussed in ( |Fang fc Wu 1996 
IRiazuelo fc Deruelle 20001 

8 (...) the general relation among the amplitudes (...) would 
be too complicated to extract any physical information out of it 
IKodama fc Sasaki 1984H . 



log T £ 
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Fig. 5. The density spectrum transfer function for per- 
turbations composed of pure growing modes. 

I24f) is an imprint of the acoustic (travelling wave) char- 
acter of the density perturbations in the radiational era. 
Consequently, the radiational era affects the structure for- 
mation processes occurring latter on, when matter domi- 
nates. Instability of perturbations is the genuine feature of 
the entire cosmological evolution rather than the property 
of separate cosmological epochs. 

6. Observables and observational constraints 

Although our primary goal is to discuss the role of de- 
caying solution in the formation of the small-scale struc- 
tures, it is worth to check the model consistence with data 
presently available for larger scales. 

Data coming from the large galaxy surveys (SDSS, 
2dFGRS and others) probe the large-scale mass distri- 
bution at late epochs, while the WMAP experiment re- 
ports the density fluctuations at the last scattering sur- 
face. These two categories of data potentially allow recon- 
structing the spectrum transfer function i.e. the scale of 
inhomogeneity enhancement after decoupling. 

To examine the spectrum transfer function (|51|l we fo- 
cus on the LSS spectral estimations, with the scale range 
of some 100 Mpc, that at the same time overlaps the 
right end of the range (k ~ few lQ^ 1 Mpc^ 1 ) probed 
by WMAP experiment. As an example we take two close 
length-scales, k ~ 0.1 Mpc and k ~ 0.05 Mpc^ 1 , for 
which the values of 8 contrast (Fig. 3 in IWu et al. 19991 
are respectively, ~ 10~ 2 and ~ 10 -3 . These scales cor- 
respond to I ~ 1800, 700 in the CMB spectrum and may 
be considered as the intermediate length scales. The spec- 
trum transfer function T t in this scales is given in the Fig. 
|H| The curves are drawn for the constant ry-surfaces be- 
tween n ~ 2 (the lowest) and 77 — 31 — corresponding 
roughly to redshift ~ 0. The increase of the model spec- 
trum transfer function l|51|l are roughly of the order of 
10 3 and 10 2 at the scales under discussion, and repro- 
duce the current density contrast, from the initial density 
contrast of the order ~ 10~ 5 . In spite of the high idealiza- 
tion, we obtain the rough agreement with observations in 
k < lMpc^ 1 regime. 
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Fig. 6. The spectrum transfer function T e as a function of 
k for diffrent rj G (2, 31) 

Substantially more difficult is to view the spectrum 
transfer function in context of extremely low, or extremely 
large scales. In the large scale limit the spectrum trans- 
fer function (|51(l becomes insensitive to changes in the 
equation of state. This confirms that the large scale per- 
turbations are "fossils" of the remnant past, what is in- 
teresting in context of their quantum or semiclassical 
origin. Recently reported (e.g. IQliveira-Costa et al. 2 003 
Hannestad 2003) discrepancy between the WMAP obser- 
vational spectrum and that coming from the inflationary 
paradigm for scales k < 0.001 Mpc -1 challenges some 
revision of the initial perturbations spectrum theory. On 
the other hand the interpretation of measurement also de- 
serves careful examination. The perturbations and the re- 
sulting temperature fluctuations are commonly related to 
each other by means of Sachs- Wolfe formula - an integral 
over gravitational potential performed along the photon 
path USachs fo Wolfe 1 9671 IPeehles 19800 . The formula is 
true in absence or neglect of the decaying modes at the 
last scattering surface IjSachs fc Wolfe 1967|) . For the den- 
sity perturbation matched to the acoustic field the decay- 
ing modes cannot be neglected, particularly not the low 
multipoles, where discrepancies become substantial 9 . To 
understand properly the measured Ci coefficients at low-1 
limit we need the exact not approximate formula for the 
temperature fluctuation. (In the discussed model the ex- 
act formula can be found, yet it is out of the scope of this 
paper.) 

Violent formation of fine scale structures in the tran- 
sient epoch, what is the characteristic feature of the model, 
enables inhomogeneities to enter nonlinear regime soon 
after decoupling. The model qualitatively supports obser- 
vations of highly developed structures at high redshifts 
z. However, WMAP angular resolution limit rules out 
any quantitative estimations for the structures less than 
10 4 M ga i. Expecting more relevant data from the forth- 

9 To neglect decaying modes one needs (&2) 2 << («2) 2 . 
After employing 13UH this condition reads 1 + 40/(3fc 2 77£ 2 ) + 
100/(fc 4 r; s 4 ) << 1, and is false for any k. Particularly, in the 
k — > the error become arbitrarily large. 



coming Planck mission we anticipate, on strength of the 
spectrum transfer function (|5I|I . that the CMBR temper- 
ature fluctuations in the fine-scale regime may have rela- 
tively low amplitude. 



7. Summary 

Generic density perturbations in the radiation dominated 
universe propagate in the same way as sounds propagate 
in air or electromagnetic waves in vacuum. As shown by 
Sachs and Wolfe ISachs fc Wolfe 19671 1 White 197311 per- 
turbations form waves travelling with the same speed v = 
independently of their scales or profiles, hence grav- 
itationally bound structures cannot form. Perturbations 
do not self-gravitate in the linear regime, so gravity may 
affect their evolution merely by affecting the dynamics 
of the homogeneous background. The wave character of 
the density perturbations (independent of their scales) 
is confirmed in the Hamiltonian formalism IjLukash 1 980 
IChibisov fc Mukhanov 1982)> . 

The hyperbolic type of propagation equation requires 
appropriate perturbation statistics, where acoustic waves 
travelling in different directions are statistically indepen- 
dent. This kind of statistics form a classical limit for quan- 
tum theories l|Lukash 19801 IChibisov fc Mukhanov 19821 
IGrishchuk 1995|l and is compatible with the gravita- 
tional waves theory fAbbott fe Harari 19861 lAllen 19961 
lAllen fc Romano 19991 |Maggiore 2"000| ). Probability and 
appropriate expectation values may depend on the wave 
frequency, but not on the direction of propagation, nei- 
ther the wave phase at any time or position. A random 
choice of plane waves guarantees that the perturbations 
and their canonical momenta are statistically indepen- 
dent and uncorrelated quantities at any time. This fi- 
nally results in the stability of homogeneous expanding 
environment. Perturbations larger than the sound hori- 
zon decay during the radiation era, while those well in- 
side the horizon keep their magnitude constant in time. 
This property, although contradicting Jeans conjecture, 
confirms results obtained in other gauges: synchronous 
HGrishchuk 19941 IGolda fc Woszczyna 200ip , longitudinal 
(Mukhanov et al. 1992 part III) and in Hamiltonian for- 
malism (Chi bisov fc Mukhanov 1982J) . Similar decrease in 
amplitude occurs for the large scale component of the 
scalar field (Stebbins & Veerarghavan 1993). 

Cosmic structure formation naturally occurs at the 
transition from the radiational to matter domination era. 
To match the acoustic field at the decoupling, the grow- 
ing and decaying modes contribute in the short wave limit 
with nearly opposite phases, and therefore, compensate 
each other at E. After the transition both modes "de- 
couple" and the resulting superposition grows explosively. 
Physically it means that the pressure decay from p = e/3 
to p — excite perturbations much lower than hori- 
zon scale, while leave untouched the amplitude of those, 
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which substantially exceed horizon 10 . For growth of inho- 
mogeneities the break down of the acoustic approxima- 
tion is responsible, i.e. the same group of physical phe- 
nomena that may excite shock waves in the interstellar 
medium. This kind of instabilities cannot be described by 
any formalism, which a' priori neglects the role of de- 
caying modes, no matter how realistic are the models for 
physics of recombination and decoupling, which are used. 

This is obvious that the sharp transition between cos- 
mological eras is not a realistic model for decoupling or re- 
combination. In reality these phenomena are continuous, 
take some cosmologically substantial time, and involve a 
number of complex physical processes. The recombination 
and decoupling do not coincide. One can hardly expect 
that the realistic situation can be described by simple an- 
alytical solutions as presented in this paper. Eventually 
one have to apply numerical codes involving multi fluid 
hydrodynamics or magnetohydodynamics. Still, the prob- 
lem of the initial state remains. The physical meaning of 
the obtained numerical results strongly depends on their 
stability against initial condition and on the physical rel- 
evance of the initial state assumed at early epochs. Most 
of hydrodynamical codes are ready to work with travel- 
ling waves, therefore, the complete numerical analysis of 
random acoustic fields in the expanding universe- without 
neglecting a' priori the role of 'decaying modes' - is ba- 
sically possible. On the other hand, simple but nontrivial 
analytic solutions presented in this paper may easily be 
used to verify numeric procedures. 
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